Construction of a giant vortex state in a trapped Fermi system 
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A superfluid atomic Fermi system may support a giant vortex if the trapping potential is anhar- 
monic. In such a potential, the single-particle spectrum has a positive curvature as a function of 
angular momentum. A tractable model is put up in which the lowest and next lowest Landau levels 
are occupied. Different parameter regimes are identified and characterized. Due to the dependence 
of the interaction on angular momentum quantum number, the Cooper pairing is at its strongest 
not only close to the Fermi level, but also close to the energy minimum. It is shown that the gas 
is superfluid in the interior of the toroidal density distribution and normal in the outer regions. 
Furthermore, the pairing may give rise to a localized density depression in configuration space. 
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■ I- INTRODUCTION 



5— i ' Among the features that distinguish a superfluid from a normal fluid is the quantization of fluid circulation. This 
implies that a superfluid responds to external rotation by forming one or more vortices, each of which carries an integer 
number of circulation quanta. In an infinite homogeneous fluid, it is well known that singly quantized vortices are 
energetically favorable for a given rotation rate^since the angular momentum is proportional to the quantum number 
and the energy is proportional to its square [HQ- As a result, a rotating homogeneous superfluid does normally not 
form giant vortices, i. e., vortices with quantum number larger than one, but instead a lattice of singly quantized 
vortices will form. Such lattices are also observed in bosonic H H and fermionic [I| superfluid gases contained in 
harmonic trapping potentials. The parameters in these experiments are such that the system is much larger than the 
core of a single vortex, so that the energy minimization argument sketched above goes through. 

The situation is different if the length scale of a vortex is comparable to the system radius. In finite superconductors, 
it has been found numerically that giant vortices form in small samples at strong fields §]. For a Bose condensate, 
the quantum number of a vortex depends on the shape of the potential used to confine the atoms as well as on 
I ' parameters associated with rotation and inter-particle interactions. When the potential in the radial direction is 
^ . harmonic, an array of singly quantized vortices is the energetically favorable state for all parameter values [HH, but 
CSJ ' in an anharmonic potential, a Bose condensate will dev elop a gian t vortex if the rotation speed is high enough and 
the interaction energy is weak enough 0, El El 13 13 El El 111 III E3- 

Because a condensed Bose gas at zero temperature obeys a nonlinear Schrodinger equation that in the noninteracting 
limit reduces to the single-particle Schrodinger equation [2, H| , it is possible to put up an exact argument proving the 
■ energetic stability of giant vortices in anharmonic traps and their instability in harmonic traps In a Fermi gas, 
no such criteria have yet been established. Therefore, as a first step towards understanding the rotational properties 
. in superfluid Fermi gases, this paper shows that it is possible, within a BCS ansatz, to construct a situation where 
a giant vortex is the energetically favorable state. In order to make the problem tractable, the study is confined to 
the so-called lowest Landau level (LLL), i. e., the space of states without radial excitations; the extension to the next 
lowest Landau level will also be studied. In Sec. [n] the physical setting is described and definitions are introduced. 
In Sec. IIHI the solution of the equations is carried out. Section llVl discusses the physics of the solution in different 
parameter regimes. In Sec. \V\ the effects of extending the space of occupied states and approaching more realistic 
parameter values are discussed. Finally, Sec. IVII provides a conclusion and outlook. 



II. LOWEST LANDAU LEVEL 

Consider a system of spin-1/2 fermions confined in an anharmonic, i. e. quadratic plus quartic, trapping potential. 
The fermions are subject to a rotational force with a rotation frequency fl. (Alternatively, Q can be seen as a Lagrange 
multiplier applied to ensure a finite angular momentum L z ). The number of particles in each spin state is N, so that 
the total number is 2N. The confinement in the direction along the rotation axis is assumed to be so tight that the 
motion in that direction is frozen out, so that the system can be assumed to be two-dimensional. The Hamiltonian is 

H = J2 I dnl>t(r)H Mr) + 9 [ ^|(r)^(r)^(r)V T (r), (1) 



2 




FIG. 1: Energy levels in the two lowest Landau levels in an anharmonic trap with external rotation frequency SI = 1.1 
and anharmonicity a — 0.01. Squares are the exact numerically computed energy levels. Asterisks denote the variational 
approximation using harmonic-oscillator eigenfunctions. 



where a runs over the two spin indices 1,1, and the interaction strength g is assumed negative, g = 
units so that the particle mass and Planck's constant Ti are unity, the single-particle Hamiltonian is 

H = -lv 2 + v(r)~nL z . 



\g\. Choosing 



(2) 



The potential is 



V(r) = — r 2 (l + ar 2 ), 



(3) 



so that u> is the trap frequency associated with the harmonic part and a is the degree of anharmonicity. The passage to 
dimensionless units is completed by setting w = 1. The single-particle spectrum associated with the Hamiltonian (J2J 
can be solved numerically, resulting in a positive curvature for the energy as a function of angular quantum number m 
for a fixed radial quantum number n r , as seen in Fig.^ In order to make the problem tractable one needs to assume 
that no states are occupied except those in the lowest Landau level (LLL), i. e. the levels with no radial nodes, n r = 0. 
The restrictions on physical parameters needed for this assumption to hold will shortly be explored. In general, the 
single-particle spectrum in the LLL around the minimum can be approximated with a quadratic expression, 



(4) 
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where k = to — mo is a shifted azimuthal quantum number; the offset Too is the quantum number at which the energy 
has its minimum; and the curvature a has to, in generai, be computed numerically. When the anharmonicity a is 
very small, the energy spectrum can be solved perturbatively using harmonic-oscillator eigenfunctions, resulting in 

e m = ui - (Q- l)m+ ~(m+l)(m + 2). (5) 

Thus, in this limit the minimum-energy quantum number is given by 

m = (ft- l)/a-3/2, (6) 

the curvature is a — a, and the ground-state energy eo = 1 — to,q + a. For not so small a, or high enough quantum 
numbers to, the corrections to the perturbative expressions become important, but they only result in quantitative 
shifts. In the following, we shall use the effective curvature parameter a, which may be slightly different from a, in 
the calculations, but the more physical discussions and order-of-magnitude estimates will be stated in terms of the 
actual anharmonicity a. 

The number of particles is determined by the chemical potential /i, which in the weakly interacting limit is equal 
to the Fermi energy. For convenience, a Fermi quantum number up is defined through the relation 

A* = e o + ( 7 ) 

so that the actual number of particles per spin state is N = 2kf in the weakly interacting limit, but one must bear 
in mind that it may differ from that when interactions are finite. The parameters are assumed to be such that /i < 0, 
implying too — up > 0. That way, the occupied part of the spectrum is bounded from below at some finite angular 
momentum quantum number m > 0, and one can treat the spectrum as if it is extended indefinitely in both the 
positive and negative K-direction. 

The coupling constant q introduced in Eq. is a pseudopotential whose relation to the actual scattering length a s 
is in general nontrivial 20], but in the weak-coupling limit it reduces to the simple form g = Ana s . The conventional 
dimensionless coupling parameter in the study of trapped Fermi gases is the product Q = kpa s , the Fermi wavenumber 
times the s-wave scattering length. In terms of our Fermi quantum number Kp and coupling constant g, we obtain 

= a 1 / 2 Ki?g/47r, and we shall confine the study to small values of £, i. c. weak coupling, in order not to violate the 
conditions for the model to be valid. We now discuss the conditions for the LLL assumption to hold. Firstly, the 
next lowest Landau level, with n r — 1 , must lie higher than the chemical potential. The spacing between the Landau 
levels is 2 (in the present units where u> = 1), so the condition is /i < eo + 2 or in terms of physical quantitites, 

a {N\ 2 „ 

Furthermore the condition [i < stated above needs to be improved. Interactions broaden the occupation of fcrmions 
and smear out the Fermi level on a scale |g|, the interaction strength. The distance between the lower Fermi level and 
the edge of the spectrum at to = thus needs to be larger than \g\, that is, to — up > \g\ or 

Q — 1 AT 

T >5- (9) 

a I 

Note that according to the first inequality, Eq. JSJ), a must be of order N~ 2 or smaller which means that the second 
restriction © is automatically fulfilled as long as SI — 1 is positive and not smaller than \g\N~ 2 , which is an extremely 
small quantity. The first inequality does, however, put quite severe limitations on the experimental implementation 
of the state, since it forces the anharmonicity to be less than or of the order of the inverse square of the number of 
particles. These matters will be discussed further in the concluding section. 



III. BOGOLIUBOV-DE GENNES ANALYSIS 



Through a number of assumptions on the physical parameters we have arrived at a description of a Fermi gas 
with a quadratic single-particle spectrum. This is the textbook case, and in principle, one may simply apply the 
usual textbook BCS expressions known for superconductors Q. Two things remain to be sorted out: first, the gap 
must be assumed to be larger than the level spacing, which permits us to treat the levels as a continuum (for paired 
Fermi systems with a discrete spectrum, see Ref. 21]). Second, one must consider the variation of the coupling with 



4 



quantum number k. This will have an effect on the density profile and quasi-particle energy spectrum, as we shall 
see. 

The Bogoliubov-de Gennes (BdG) equation in two dimensions reads 

ffo-C AM) \(u K (r,0)\_ fu K (r,9)\ 

A(r, 9)* -HS + M J\ v K (r, 0) )~ K \ v K (r, 0) ) ' [W) 

where we have anticipated that the solutions u K , v K can be labeled by the shifted angular momentum quantum number 
k. When the single-particle energy spectrum is symmetric around m = too (or k = 0), the gap function can be written 

A(r,0) = A{r)e' i2m ° e . (ff) 

This amounts to assuming a cylindrically symmetric gap function and corresponds to the usual assumption of a 
homogeneous gap function in an infinite system. In order to match the 9 dependence, the amplitudes take on the 
form 

u K (r) = u K (j) mo+K (r), 

v K (r) = v K 4>* mo - K (r), (12) 

where <f> m are single-particle eigenfunctions in the anharmonic trap and u K and v K are space-independent amplitudes. 
The BdG equation l|l(J|l is now multiplied through by eigenfunctions <p m and integrated, resulting in 

U-- 5 ;)(:H(:;)' ™ 

where 



A K = J dr2nrA(r)^ no+K (r)cj ) * mo _ K (r), (14) 



and 



£k = £k - M = ^( k2 - k f)- ( 15 ) 



The solution for the energies and amplitudes is similar to the traditional BCS expressions: 



E K = + \A K \*, 



, 2 (l + ^-),vi = -(!-£). (16) 

Now turn to the self-consistency equation for the gap function, 

A(r) = | ff |5>«(rK(r). (17) 

K 

Inserting the calculated Bogoliubov amplitudes into the self-consistency equation, one obtains 

A(r) = | 5 |^0 mo+K (r)0 mn _ K (r)^. (18) 

Multiply by <p mo+K i (r)(j) ma - K i(r), integrate and rename the subscripts. The result is 

A « = ME§7V as) 

where the coupling matrix element was defined as 

V KtK > = (mo + K, to - k|to + k', rn - n'), and 
(m,n\k,l) = / drfafinfatf ■ (20) 
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In order to be able to produce a definite expression, let us make use of the results for a harmonic trap. Since Eq. JSJ) 
confines the parameters to lie in the regime of a very weak anharmonicity, the anharmonic trap can at most bring 
in small quantitative corrections. The matrix element for the contact potential in the space of harmonic-oscillator 
eigenf unctions is [j| 



{ m,n\k,l) = 1 [(H + M + |fej + W)/2]! 

27r2(M+M+l fc l+l i l)/VH!H!|fc|^|! 



Note that the potential is separable: 



K,k< - V K V K ,, V K y 27r22mo(m ( o + °) !(mo _ k)V ( 22 ) 

and furthermore, the potential is in the limit of large mo approximated by 

V K w — 1 -e^. (23) 

A separable potential presents a considerable simplification. The solution for the gap is 

A K = AV K , (24) 
where the constant A is found by solving the self-consistency equation in its final form, 

l = \g\Y-= % (25) 

This sum has, in general, to be computed numerically. Nevertheless, there is plenty of information to be extracted 
from its qualitative form, as we shall explore in the next section. 

IV. PARAMETER REGIMES 

Consider the self-consistency equation l(23|) . If the coupling is weak, A is expected to be small as well. In this case 
the largest contribution to the self-consistency sum comes from the states close to the Fermi surface, and the usual 
textbook BCS solution is recovered 0, 

A = w c e-VW°W»W), (26) 

where uj c is the cut-off frequency for the sum and N(0) is the density of states at the Fermi level. However, that solution 
is only valid when the interaction strength at the Fermi surface is larger than the separation of the energy levels; if 
it is smaller, one cannot neglect discreteness effects j^. This criterion translates to |g|e~ K ^/ m °/y / 7TTOo > anp- On 
the other hand, if \g\ is too large, pairing will happen away from the Fermi surface, and as will shortly be shown, this 
occurs when \g\ > anp. Combining the two inequalities, one finds that a criterion for avoiding both discreteness and 
off-Fermi surface pairing is 



>e K "l m \ (27) 



which cannot be fulfilled for any value of kf- We conclude that the standard textbook solution is not applicable in 
the present system. 

We now derive the criterion already used above, for states away from the Fermi surface to contribute to the 
self-consistency sum. When the sum is exhausted by terms far from the Fermi surface we can write 

\g\ Y . J K =1. (28) 



Setting A = yields an inequality, 



\k\<kf ^ K 



l ii >1. (29) 
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FIG. 2: Level curves for the gap function A in three slices of parameter space, (a) Fixed anharmonicity a — 0.0001, (b) fixed 
Fermi quantum number kf = 60, (c) fixed coupling \g\ = 1.5. The minimum-energy quantum number is fixed to mo = 100 in 
all three panels. Asterisks denote the phase-space points used in Fig. [3 



The sum is estimated by inserting the maximum value of the summand and multiplying it by the effective range of 
the potential V K , 



|« ? | N /2^-°- > 1, (30) 
l?o| 

where \J2m^ is the width of the function V K and V 2 = (47r 3 mo) _1/ ' 2 is its value in the region where the summand is at 
its largest. The factors involving mo cancel out, so the critical coupling strength for an off-Fermi surface contribution 
reduces to 

Iffl > 9c oc |f | oc olk%. (31) 

In Fig. |21 the numerically computed value of A is plotted in various slices of parameter space. It is seen how the 
A = 0.01 level curve can be reasonably well fitted to Eq. H31[) in all three plots. Beyond this curve, the discreteness 
of the spectrum is expected to be important and the BCS solution cannot be trusted to be quantitatively correct. 

The quasi-particle spectrum differs qualitatively between the different parameter regimes. This is illustrated in 
Fig. [3J When \g\ < g c , one recognizes the textbook situation: pairing takes place mostly at the Fermi surface and 
the spectrum is equal to the absolute value of the single-particle spectrum except close to the Fermi surface, where 
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FIG. 3: (a),(c): Full line: Quasi-particle energies E K . Dashed: K-dependent gap function A K . Dotted: summand S K in the 
self-consistency sum. The gap function is magnified by a factor 10 s in (a). Horizontal line: chemical potential fi. (b), (d): 
Spatial distributions of density (full line) and gap function A(r) (dashed). The gap function is scaled to have a peak value of 
0.4. Parameters are chosen as mo = 100, kf = 60, and a = 0.0001. (a-b): \g\ = 0.3, (c-d): \g\ = 3.0. 



a gap opens. On the scale of Fig. [3^, the gap is too small to be visible. The plot also includes the /t-dependent gap 
function A K and the summand in the self-consistency equation, Eq. (|25[) . S K — V^/E K . It is seen that an appreciable 
contribution to the sum comes from terms close to the Fermi surface. If \g\ were even smaller, the sum would be 
entirely exhausted by the terms close to the Fermi surface, but as was found above, the BCS solution would be altered 
by effects associated with the discreteness of the spectrum. For stronger coupling, \g\ > g c , the pairing happens 
predominantly close to k = 0, where the interaction is at its maximum, as seen in Figs. |3t-d. The quasi-particle 
spectrum is affected in this region, but at larger ft it resembles the usual BCS spectrum. 

The spatial profile of the gap function is completely insensitive to the profile in K-space of the summand as long as 
the population is confined to the LLL. It is given by the self-consistency equation in configuration space, Eq. (|18f) . 



Thus, the gap function is a constant times the 2mo'th radial harmonic-oscillator cigenfunction. At the same time, 
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the density is much broader as a function of r, and it does depend on the coupling. The density is 

n(r)-2^|0 mo _ K (r)| 2 K| 2 , (33) 

K 

with \vk,\ 2 = (1 — £ K /E K )/2. It is a peculiar feature of the LLL states that each eigenfunction mo _ re (r) is confined 
to a very narrow, shell-like region in configuration space, centered at r — (mo — ft) 1 / 2 . As a result, the distribution 
of particles in K-space is mirrored by the spatial density profile, giving direct observational access to the distribution 
of particles in angular momentum space. For weak coupling, all states below the Fermi surface contribute, just as in 
ordinary BCS theory. The spatial profile as a function of the radial coordinate is a flattened, almost square profile, 
nicely mirroring the filled Fermi sea. As a result, only the interior of the toroidal density distribution is supcrfluid, 
while the outer regions are normal. For stronger coupling, the Bogoliubov amplitude v K is depleted around k — and 
the density profile has a dip in the region where the gap function is at its largest. Mathematically, it is easy to see 
why the density is depressed at t he trap b ottom when the coupling is strong: the Bogoliubov amplitude v K is equal to 
half its peak value when E K = y A 2 + £ K ^ £«• Physically, it is the strong interactions that alter the occupation of 
the single-particle levels. As a result, one obtains a spatial region of strong Cooper pairing that shows up as a density 
depression. In the next section, we shall examine the conditions for this effect to persist when higher Landau levels 
are included. 



V. NLLL 



In order to approach slightly more realistic territory in parameter space, let us now study what happens when 
the next lowest Landau level (NLLL) is brought into the calculation. Inspection of Eq. (JSJ reveals that this means 
allowing the product aA 2 /16 to be slightly larger than unity. This system is still not easily realizable in experiment, 
but it is a first correction to the even more difficult LLL case studied above. The new Landau level introduces a 
new parameter into the system, namely, the energy difference between Landau levels, which is equal to 2cu — 2 in 
dimensionless units. Thus, the balance between three main parameters determine the physics of the system as well 
as the validity of the truncation of the basis. The three are the chemical potential /i, the coupling strength g, and 
the energy difference 2o>. In addition, the balance between the Fermi level kf and the width of the coupling, y/mo, 
plays a role in determining the physics. 

The single-particle eigenfunctions in the NLLL are 



</w(r) = ( 1 - ^— ) (re ie re- r ^ 2 . (34) 

V nm\ \ m + 1 J 

The perturbative expression for the energies in the quartic trap are (see Fig. ^| 

e m ,i = 3oj - (Cl - Lo)m + -(m + 2)(m + 4). (35) 

The position of the energy minimum of the NLLL coincides with that of the LLL, mo, to within 3/2 quanta. That 
shift can safely be neglected here. On the other hand, when many Landau levels are brought into the calculation, the 
resulting asymmetry of the single-particle energy spectrum is likely to break the cylindrical symmetry such that the 
order parameter becomes a superposition of several angular-momentum states. The resulting profile is expected to 
be an annulus pierced by a ring of singly quantized vortices [l4| . 
Now consider the coupling matrix elements. They are 

^0001 



VZ = ((mo + K,0),(mo-K,0)\(m o + q,0),(mo-q,l)) 
= (2m )!(l-2g) 

27r2 2m " ^(m Q + q+ l)!(m -q + l)!(m + n)\(m - k)\ 



( K + q )/(2mo) + 0( m -2) ! ( 36 ) 



TfOOll 

K,q 



1 (2m )!(m + 1 - 2q 2 ) 

47T 2 2m ° y/(m + q+ l)!(m -q + l)!(m + n)\(m - k)1 ' 
f-2 g 2 /m c _ (K 2 +t> 2 )/(2mo ^ 

47Ty / 7rmo 



(37) 
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(2m )!(l + m 



Anq) 



87r2 2m o y/(m Q + q)\(m -q + l)!(m + «)!(m 

1 _^2x„2^/rt m „l A 4k<? 



1)! 



-(» +9 )/(2m ) x + 



r 1 



1101 



47Ty / 7rmo ^ 2mo 

(2mo)!(l + mo — 2q — 2m q — 2k 2 



AK 2 q) 



2?r2 2m « y(mo + g )!(m - 

1 ~ ^_ p -(K 2 +g 2 )/(2m ) 

rn 



l)!(m 



l)!(m 



1)! 



2tt 3 / 2 

1 (2to )![2 + 3m 2 , - 2(k 2 + q 2 ) + to (5 - 2(« 2 + q 2 )) + 4«V] 
8ir 2 2m ° y/(m + q+ l)!(m - g + 1)!(to + n + l)!(m - k + 1)! 



87Ty / ¥mo 



(1- 



l + 2( K 2 + g 2 ) )c „ (K 2 +q2)/(2mo) 



3mo 



(38) 



(39) 



(40) 



The other matrix elements can be obtained from these by a sign change, e. g. V®° q 10 = V® 00 ^. It is seen that these 
new matrix elements are not separable in the way that simplified the calculations within the LLL. Moreover, the 
matrix elements I/ 0001 and V 1101 are smaller than the others by a factor mj 1 / 2 ; however, that does not mean that the 
inter-level coupling cannot be neglected, since the matrix elements V 0011 and V 0101 are still as large as the intra-level 
matrix elements. 

Employing the cylindrically symmetric assumption for A(r) (keeping in mind that this may change when more 
Landau levels are brought into the problem), the ansatz for the Bogoliubov amplitudes has to be labeled by an 
angular momentum quantum number n and a branch index j = 0, 1 (enumerating the two solutions that exist for 
every k), 



v Kj (r) = v Kj0 (t>* mo - K ,o( r ) 
The Bogoliubov equation for the (k, j)'th mode is 



MKjl0m o + K ,l( r ) ; 

-VKj\4>* mo - K ,i{r). 



( ?k0 



A 00 * 

v Air 



A°° 

Ar -u 

A"* 



A" 1 
A" 






( U «J0 \ 




( u *j0 \ 


) 


U K jl 


= E Kj 


U K jl 

v Kj0 1 




V V K jl ) 




V V K ji ) 



(41) 



(42) 



We have defined £ K 



H and the matrix elements 



A"™ = / d 2 rA( 



The self-consistency equation is in general 



and integrating one obtains 



A(r) - \g\J2u Kj (r)v Kj {r)*, 



\g\ U J«'«3 W j* K 'n 4 

K f j,Tl3 ,714 



(43) 



(44) 



(45) 



The coupled equations are solved numerically. Figure ^displays a few illustrative cases. It has been checked in trial 
calculations that the inclusion of more Landau levels does not alter the density profiles shown in this paper, and only 
marginally distorts the quasi-particle energy curves. Basically, there are four parameter regimes of interest. When 
fi <C 1 and \g\ <C 1, the NLLL is not occupied at all. This case was studied in the previous section. When fj, ~ 1 
and \g\ <C 1, the Fermi surface cuts through both energy bands, but the coupling is not too strong, as in Figs.^-b. 
There is a gap at the bottom of each of the two bands. The NLLL has a larger gap because the coupling is stronger 
for smaller k. The associated density profile is expected to be bimodal, with an enhanced density in the interior 
of the toroidal cloud where two Landau levels are populated (cf. [13 )• Because the Fermi level has been limited to 
Kp = 60 in the numerical calculation, this "wedding cake" structure does not come out as clearly in Fig. ^3 as it 
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-60 -40 -20 20 40 60 5 10 15 20 

k r 



FIG. 4: Energy spectra and density profiles calculated using two Landau levels (NLLL approximation). (a),(c),(e): Quasi- 
particle energies E K . Horizontal line: chemical potential fi. (b), (d), (f): Spatial distributions of density (full line) and gap 
function A(r) (dashed). The gap function is scaled to have a peak value of 0.4 times the maximum density. Parameters are 
chosen as m = 100, and n F = 60. (a-b): \g\ = 0.4 and a = 0.0012, (c-d): \g\ = 8.0 and a = 0.0001, (e-f): \g\ = 8.0 and 
a = 0.001. 
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would in a bigger system. Figs. 0J>d illustrates the case where the chemical potential is small, but the coupling is 
comparable to the separation between the energy bands; fj, <C 1 and \g\ ~ 1. The Fermi level lies far below the NLLL, 
but because of the strong coupling, both levels are deformed. The resulting density profile is still depleted in the 
region of strong Cooper pairing, because the contribution from the NLLL is modest. For stronger interactions, more 
Landau levels are coupled and the NLLL approximation would fail. Note that the effect of the strong coupling on 
the quasi-particle spectrum is much more pronounced than the effect on the density. Finally, in Fig.^Jj-f, both p, and 
\g\ are comparable to the inter-level gap. The energy levels are heavily distorted by the strong coupling. The density 

1/2 

near r = m Q is no longer depleted, because the NLLL states now occupy that regime. The gap function takes on 
a slightly nontrivial shape in configuration space; this is again an effect of the relatively few single-particle levels 
involved in the calculation. Trial calculations indicate that when the number of states is increased, the gap function 
acquires a smoother shape, but the contribution from the radially excited states broadens its profile compared to the 
LLL case. 

Extrapolating from the results presented here, one may observe that the density is affected by two competing 
effects: as the Fermi surface includes more and more Landau levels, a "wedding cake" structure is expected in the 
density for weak interactions, but at the same time, strong interactions will deplete the interior of the toroidal density 
distribution. 



VI. CONCLUSION AND OUTLOOK 



This paper has shown that a Fermi superfluid trapped in an anharmonic potential can develop a giant vortex if it is 
rotated at high enough angular velocity, since the single-particle energy spectrum will then have its global minimum 
at a finite angular momentum. A giant-vortex state has been constructed in the basis of the lowest radially excited 
states, the lowest Landau level (LLL). Although attaining the LLL regime experimentally puts severe constraints on 
precision, the theoretical construction of such a state demonstrates that there exist parameter regimes in anharmonic 
traps for which a Fermi gas develops a giant vortex. As long as only the LLL is occupied, the density and gap function 
are cylindrically symmetric and form a toroidal profile. The superfluid region is confined to the interior of the torus, 
while the edges are normal. Since Cooper pairing is strongest in the center of the toroidal profile, the density is 
depleted in that region. The effect on the density appears, however, to be exclusive to the LLL. As more Landau 
levels are occupied, singly quantized vortices may form within the torus. 

In order to have only the lowest Landau level populated in the quadratic plus quartic trap, the anharmonicity and 
rotation frequency have to be extremely fine tuned. In addition, any mechanical stirring method Q will bring in 
corrections to the anharmonic trap that may be hard to control. It is therefore probably better to control the angular 
momentum of the gas by first stirring it in a harmonic potential and subsequently applying the anharmonicity. As an 
alternative to an extremely small quartic trapping term, an optical hard-wall potential at a suitably large distance 
from the trap center may be a more practical choice. It is possible that such a setup, with extensive amounts of fine 
tuning, may be a more practical route towards accessing the lowest Landau level. 

Nevertheless, this study should first and foremost be seen as a demonstration by construction that trapped Fermi 
superfluids can sustain giant vortices. Several important issues remain to be studied, notably the size and shape of 
the allowed parameter regime for giant vortices; the physics in the BEC-BCS crossover regime; and whether giant 
vortices can ever be the favorable rotational configuration in purely harmonic traps. We hope to address these issues 
in future work. 
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